The establishment of ecological conservation for herpetofauna species in hotspot areas of South Korea

Understanding the geographic distribution of species is crucial for establishing protected areas. This study aimed to identify the preferred habitat environment of South Korean herpetofauna using distribution point information, providing the information necessary to protect their habitat by establishing a species distribution model. We found that climate variables in the region where 19 amphibians and 20 reptiles were distributed correlated with the altitude, suggesting that altitude had a major influence on their distribution. The species distribution modeling indicated that 10–12 amphibian and 13–16 reptile species inhabit the Gangwon-do region, forming hotspot areas in the eastern and western regions around the Taebaek Mountains. Some of these hotspot areas occurred in the Demilitarized Zone and national parks, which are government-managed ecological conservation areas. However, some hotspot areas are vulnerable to habitat destruction due to development and deforestation as they are not designated conservation areas. Therefore, it is necessary to establish new conservation areas with a focus on herpetofauna after confirming the actual inhabitation of species through precise monitoring in predicted hotspot areas and designating them as protected areas. Our results can serve as important basic data for establishing protection measures and designating protected areas for herpetofauna species.


Results
Habitat distribution characteristics. The 19 amphibian and 20 reptile species inhabiting South Korea were observed at 25,400 and 8,581 locations, respectively (Table 1). Rana nigromaculata was the most commonly identified (5306 locations), and Hynobius yangi was the least commonly identified (21 locations) amphibian species, and Rhabdophis tigrina was the most commonly identified (2121 locations), and Eremias argus was the least commonly identified (212 locations) reptile species (Table 1).
Species distribution model. The average AUC of the 19 amphibian and 20 reptile species was 0.798 ± 0.139 and 0.764 ± 0.103, respectively. The 10% training error value of the test sample was 0.120 ± 0.020 for amphibians and 0.180 ± 0.074 for reptiles, indicating that the overall error was low ( Table 2).
The variable with the highest contribution in the distribution model of amphibians was altitude (39.10%), followed by the variables Bio1 (16.71%) and Bio4 (14.46%) (Fig. 2). Altitude showed the highest contribution in 12 out of 19 models ( Table 2). The variable with the highest contribution in the distribution model of reptiles was also altitude (25.79%), followed by the variables Bio2 (18.67%) and Bio1 (11.03%) (Fig. 2). Altitude showed the highest contribution in 10 out of 20 models (Table 2).

Main distribution areas.
Herpetofauna species were predicted to be distributed in most areas of South Korea. For amphibians, the areas where 7-9 species coexisted were predicted to be the most with 72,449 cells, and for reptiles, the areas where 9-12 species coexisted were predicted to be the most with 59,728 cells. The hotspot areas for amphibians, where the most diverse species coexisted, accounted for 26,434 cells with 10-12 species, and the hotspot areas for reptiles accounted for 7,823 cells with 13-16 species (Figs. 3a,b, 4a,b). The hotspot areas of amphibians included Gangwon-do (23.80%), Chungcheongnam-do (22.39%), and Jeollabuk-do (11.18%), and the core distribution areas of reptiles included Gangwon-do (29.55%), Gyeonggi-do (26.58%), and Gyeongsangnam-do (17.18%) (Figs. 3a,b, 4a,b).
In the 10,169 cells designated as national conservation areas, the areas where 7-9 amphibian species coexisted were predicted to account for the most, at 32.44%, while the areas where 9-12 reptile species coexisted were predicted to account for the most, at 31.23%. Hotspot areas in the national conservation areas accounted for 19.34% for amphibians and 6.47% for reptiles (Fig. 5).

Discussion
In the present study, major habitats of herpetofauna species were predicted using distribution models of 19 amphibian and 20 reptile species inhabiting South Korea. A high correlation was identified between important climatic variables in the areas where herpetofauna species were distributed and geographical variables, and altitude was found to be an environmental variable with the most influence on their distribution. Hotspot area predictions showed that the province with the most diverse species was Gangwon-do, around the Taebaek Mountains.
Altitude was an important environmental factor affecting the distribution of herpetofauna species, showing a high contribution in most models. In general, different herpetofauna species are distributed at different altitudes  25,30,35 . Most of the 20 reptile species are mainly distributed at altitudes between 0 and 500 m, except for the following three species: Eremias argus, which is known to prefer low altitudes between 0 and 100 m, and Gloydius intermedius and Scincella huanrenensis, which are known to prefer altitudes higher than 400 m 7,29,34 . The results of the present study were consistent with the elevation distributions for herpetofauna species reported in previous studies. We found that herpetofauna species distributed in South Korea preferred forests and farmland, with amphibians more commonly inhabiting farmland and reptiles more commonly inhabiting forests (Fig. 1b). Paddy wetlands provide an essential aquatic environment for skin-breathing amphibians, and 16 of the 19 amphibian species inhabiting South Korea, except for Kaloula borealis, Rana plancyi, and Karsenia koreana, use paddy wetlands as breeding grounds 28,30,[36][37][38] . Furthermore, among these 16 amphibian species, all but two species, Rana plancyi and Hyla suwonesis, are known to prefer paddy wetlands located near mountainous areas rather than plains, and their population size is also known to be larger than that of the other two species 27,30,38,39 . Among the 20 reptile species, all but two species, Eremias argus and Gekko japonicus, prefer forests, using rivers, valleys, ridges, grasslands, wetlands, and other habitats in the surrounding areas for breeding and hibernation 29,34,[40][41][42][43][44] . Furthermore, 18 out of 20 reptile species, excluding Gloydius intermedius and Scincella huanrenensis, appear in paddy wetlands, which they use for foraging 29,34,40 . Therefore, forests and agricultural land, which were the main habitats of herpetofauna species identified from a macroscopic point of view, adequately reflected the main habitats of herpetofauna species reported in previous studies. The most important habitats were paddy wetlands for amphibians and forests for reptiles 27,31 .
The geographic distributions of herpetofauna species predicted by species distribution modeling were consistent with the geographic ranges reported by previous surveys (Supplementary Figs. S1, S2). According to previous studies, nine amphibian species are widely observed inland, while the others have a limited distribution range 24,30,33 . Onychodactylus fischeri, Bufo stejnegeri, and Rana huanrenensis are densely distributed in the northeastern regions, including Gangwon-do, Gyeonggi-do, and Gyeongsangbuk-do. Kaloula borealis, Rana plancyi, and Hyla suwonesis are found in central and western regions, including Chungcheongnam-do and Gyeonggido, while Karsenia koreana is found in the central regions of Chungcheongnam-do and Chungcheongbuk-do. Rana catesbeiana is found in southern regions, including Gyeongsangnam-do, Jeollanam-do, and Jeju-do, and Hynobius yangi is concentrated in Gyeongsangnam-do 26,39,40,45 . Regarding reptiles, 13 species are widely observed inland, while Scincella huanrenensis and Gloydius intermedius are mainly distributed in northeastern regions such as Gangwon-do, Gyeonggi-do, and Gyeongsangbuk-do; Eremias argus in central and western regions such as Chungcheongnam-do and Gyeonggi-do; and Gekko japonicus, Sibynophis chinensis, Coluber spinalis, and Amphiesma vibakari ruthveni in southern regions such as Gyeongsangnam-do, Jeollanam-do, and Jeju-do 40,46,47 .
Recent studies have indicated the need for using habitat prediction models to establish protected areas for wild animals and plants in South Korea 5,33,48,49 . Habitat prediction model studies can be used to provide objective and scientific methods and procedures for the establishment of protected areas 4,6,20 . For example, hotspot areas were identified, and protected areas were established or proposed for 16 Hylidae species that inhabited South America, 7 Viperidae species in Africa, and all herpetofauna species in Madagascar and Morocco 4,6,20,50 . The main distribution areas of the seven endangered herpetofauna species in South Korea were found to be Chungcheong-do and western Gyeonggi-do, and the main distribution areas of three amphibian species of the genus Rana and three reptile species of the genus Gloydius are located around the Taebaek Mountains in Gangwon-do 7,29,30 . These areas were similar in location and extent to the hotspot areas identified in the present study and included national conservation areas such as national parks. Nevertheless, some modifications should be made, and additional conservation areas should be established considering the hotspot areas where many herpetofauna species were observed.  www.nature.com/scientificreports/ In the present study, additional hotspot areas for herpetofauna species were predicted to be located in Gangwon-do than in the other provinces (Fig. 3). In Gangwon-do, forests account for 81% of the total area, mostly because of the presence of the Taebaek Mountains. As Gangwon-do is the administrative district with the lowest population density, various ecological conservation areas designated by the government are located here 51,52 . The national conservation areas in Gangwon-do, where many of the hotspot areas for herpetofauna species are located, include the Demilitarized Zone (DMZ) located in the northern region; the water resources conservation areas of Chuncheonho Lake, Paroho Lake, and Soyangho Lake located in the western region; and Seoraksan and Odaesan National Parks located in the eastern region (Fig. 4). In these areas, the ecosystem is not damaged, and high biodiversity is maintained as access by civilians has been restricted since the end of the Korean War in 1953 53,54 . According to previous surveys conducted in the DMZ from 1989 to 2016, 16 amphibian and 18 reptile species were known to inhabit the area 55 . The habitat status data collected from 1997 to 2019 showed 8-12 amphibian species and 10-14 reptile species in the national parks located in Gangwon-do 56 . Additionally, we identified hotspot areas not only in the conservation areas but also in other areas within Gangwon-do due to its low population density and well-conserved ecological environments.
In the present study, the hotspot areas for herpetofauna species were identified in various regions other than Gangwon-do. In particular, hotspot areas for amphibians were mostly concentrated in Chungcheongnam-do and its surrounding areas, and hotspot areas for reptiles were concentrated in the areas around the southwest coast and the areas of the four major rivers (Fig. 4). Chungcheongnam-do and its surrounding areas consist of plains with low altitudes. In this province, well-developed paddy wetlands are used by amphibians as their main  www.nature.com/scientificreports/ feeding and breeding grounds. Unlike other regions, this province has high biodiversity, with dense populations of endangered amphibians such as Kaloula borealis, Rana plancyi, and Hyla suwonesis 7,27 . In the areas around the southwest coast and the areas of the four major rivers (Hangang River, Geumgang River, Nakdonggang River, and Yeongsangang River), high biodiversity may have been identified because of the habitats of Eremias argus, Pelodiscus sinensis, Chinemys reevesii, and Trachemys scripta elegans, which mainly inhabit coastal sand dunes, rivers, and streams, which are not present in other regions 7,40,57 . However, many of the hotspot areas for amphibians were not included in the national conservation areas.
The national conservation areas of South Korea included more than 30% of areas inhabited by 7-9 amphibian species and nine to 13 reptile species together, as well as some hotspot areas for herpetofauna species, playing an important role in habitat conservation. However, more hotspot areas were identified outside the national conservation areas. Therefore, we showed that habitat protection is not carried out in these hotspot areas not included in the national conservation areas, with a high risk of habitat destruction because of development activities such as road and apartment construction. Therefore, to protect herpetofauna species, it is necessary to establish new conservation areas focusing on herpetofauna species after confirming the actual inhabitation of species through precise monitoring in the predicted hotspot areas. Furthermore, the hotspot areas where the actual habitats were  A comprehensive analysis of the hotspot areas of 39 herpetofauna species revealed that Gangwon-do was the province with the highest number of hotspot areas, predicted to be inhabited by 10-12 amphibian species and 13-16 reptile species. In addition, amphibians were predicted to be concentrated in paddy wetlands around Chungcheongnam-do, and reptiles in the areas around the southwestern coast and areas of the four major rivers. Some hotspot areas were included within the national conservation areas, but many hotspot areas were located in areas not designated as conservation areas. It is necessary to protect the habitats of herpetofauna species by expanding the conservation areas after verification through detailed surveys in these areas. In the future, if the exact distribution range of Dryophytes flaviventris, a recently discovered species not included in this study, is revealed, more hotspot areas in addition to those revealed in this study may be discovered.

Methods
Study area. The study was conducted in South Korea, covering several regions of the Korean Peninsula and several islands, including Jeju Island. About 70% of the investigated area was covered with forests, and about 30% was agricultural land. The eastern region has a high altitude because of the Taebaek Mountains, whereas the western region is characterized by low-altitude terrain with plains and arable land (Fig. 6b,c). Korea has a continental climate with four distinct seasons, with cold and dry winters and hot and humid summers. It is divided into nine provinces: Gyeonggi-do (GG), Gangwon-do (GW), Chungcheongnam-do (CN), Chungcheongbukdo (CB), Gyeongsangnam-do (GN), Gyeongsangbuk-do (GB), Jeollanam-do (JN), Jeollabuk-do (JB), and Jeju Island (JJ) (Fig. 6a).

Species distribution data and environment analysis. The observation data of herpetofauna species
inhabiting South Korea was obtained from three sources: data from the National Natural Environment Survey conducted by the National Institute of Ecology and the National Academy of Environmental Sciences (Ref. 58 Table 1). The location of species observed in most of South Korea (98.8% of the total land area, or 99,000 km 2 ) could be confirmed based on the results of the surveys described above. All applied data were collected through field surveys by herpetofauna experts with over ten years of experience. The survey period was from early spring (February) to early winter (November), when reptiles and amphibians are active in South Korea. Its geographic scope covered the entire country, including the land and many islands 40 . All experts visually identified the species of individuals detected while walking or traveling in a car and collected geographical information. A total of 19 amphibian species and 20 reptile species were used for the analysis, excluding Dryophytes flaviventris, which was recently identified to inhabit South Korea 62 .  (Table 3). All grids were of a uniform size of 30′′ (about 1 km 2 ). In order to identify the types of habitats preferred by the studied species, the land cover map was divided into four habitat types by determining similar or overlapping variables (Supplementary Table S1). Since climate variables are highly correlated with each other, the variables with high correlation (Pearson's correlation coefficients (r) > 0.8) were excluded from the analysis to minimize the effect of multicollinearity 35,65,66 . Accordingly, the following six climate variables were used in the present study: annual average temperature (Bio1), average diurnal temperature range (Bio2), isotherm (Bio3), annual average precipitation (Bio12), summer precipitation (Bio13), and winter precipitation (Bio14). The distribution points of the herpetofauna species were projected onto all environmental variables, and the habitat environment was checked using the extracted values, and a kernel probability density plot was generated for the altitude values. According to data normality, all data were expressed as means with standard deviations, medians, or first-third quartiles. Statistical analysis was performed using R version 3.0.2 67 .

Species distribution modeling.
The maximum entropy approach model (Maxent version 3.3.3 k) 68 , one of the species distribution model (SDM) algorithms, is the most widely used for wild organisms and provides the highest prediction result based on regression analysis 33,69,70 Unlike other algorithms (e.g., GLM, GAM, RF, etc.), this model integrates pseudo-absence points without any assumption of certainty and maintains the most possible uniform distribution under the limitations imposed by predictor variables, leading to the least bias for the presence of predicted results and its most conservative estimates [71][72][73][74] . Since Maxent can predict even with a small number of samples due to the use of appearance data alone, it is actively used in studies on reptiles and amphibians that are difficult for field observation 29,[75][76][77][78] . An SDM was constructed using the appearance data of reptiles and amphibians as dependent variables and environmental variables (a total of eight environmental variables; Table 3) as independent variables. The models were repeatedly run 15 times using default parameters, including logistic output, 1 for regularization multiplier, and 10,000 for background points 29,79 . This study generated a potential dichotomous (presence/absence) distribution raster based on 10% training presence (including 90% of occurrences) as the threshold [7][8][9][10] . The 10th percentile threshold has the advantages of being less sensitive to extreme environmental values and reducing commission errors 80,81 . Thus, it is used for wild animals, including  www.nature.com/scientificreports/ amphibians and reptiles that are mobile 77,78,82,83 . To evaluate the model, the dataset was divided into a training set for 75% and a testing set for 25% through the random test percentage, and subjected to 5000 iterations 29,79 . The explanatory power of the model was verified by calculating the area under the curve (AUC), which is the value of the lower area of the curve, by receiver operating characteristics (ROC) verification. AUC values range from 0.0 to 1.0, with a value closer to 1.0 indicating a higher prediction accuracy of the model 84 . Most studies using habitat prediction programs used AUC values to evaluate model performance, which can be sensitively affected by model conditions such as the number of samples and resolution [85][86][87][88][89] . In order to supplement this, some studies related to SDM suggest an omission rate in addition to the AUC value 10,30,79 . The omission rate is calculated as the ratio of points that were not predicted based on a threshold and were thus missing. The values range from 0.0 to 1.0, with a lower value indicating fewer omissions in the analysis process. Therefore, in the present study, besides the AUC value, the omission rate shown in the 10% training presence was also considered 86 . The contribution of each environmental variable to the areas where the herpetofauna species were distributed was calculated using the average percent contribution determined by the jackknife test. The derived SDMs were overlapped based on species through Arc GIS (Ver. 10.3; ESRI, Redlands, CA, USA), and the number of cells present in each of the nine administrative areas was identified by a percentage based on interspecies overlapping areas. Furthermore, by superimposing the generated species distribution map, the geographic range of the hotspot areas where herpetofauna species were densely distributed was identified. The areas predicted to be inhabited by 10-12 amphibian species and those predicted to be inhabited by 13-16 reptile species were determined as hotspot areas 90,91 . To identify the existing national conservation areas designated in South Korea, the geographic scope of the environmental conservation areas and national parks obtained from the National Geographic Information Institute was used (Fig. 6a) 52,[92][93][94] .

Data availability
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.